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TITLE 

TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 

PARAMETERS 

ACCURACY 



DURATION 



READ AROUND 



DESCRIPTION 



A search for the Real Roots of f (x) = (DOI or SADOl) 

Closed 

80 

7 to l6 inclusive 

Several Program parameters must be set. 

Depends upon the function f (x) and the error in the routine 

which computes it a Possible error inherent in this routine 

is about 2 e 

[2 + (2.33 + 2t)n + (6 4- t)ir] ms where t is the duration 

in ms of the routine which computes f (x)« n is the number 

of intervals examined for roots % i is the average number of 

iterative linear interpolations necessary to find a root, 

and r is the number of roots found. 

Depends upon routine used to compute f (x), whose duration 

is denoted by -t(iBsJ. RAR adjacent to routine is 22 « 5/ (2 » 33 + 2t) 

or 22.5/(6+t) whichever is greater Next to the temporary 

store it may be 8 times as great * 

This routine is a combination of Code HI (inverse Interpolation)^ 

Code P2 (Decimal Fraction Print )_, plus orders which examine 

an arbitrary range of the independent variable for roots of 

f (x) - 0o In n steps of length hy the routine examines the 

interval x n < x < x~ + nh for roots ° The ith interval is 



regarded as x 







[i-l.)h < x < x„-f ihc If such an interval 



contains an even number of roots o (e.g., 0) none will be 
found j if it contains an odd number of roots (e.g. ; l) one 
will be found „ This means the programmer must have sufficient 
information about f (x) to choose x ,% and h compatible with 
both economy of time and the finding of all the roots » Multiple 



-2- 



roots are treated as single roots. 

*he routine is entered with x Q in 1^ and the orders* 

JO KF or 50 rF 
50 pF 



p+1 



p+2 



p+5 



26 qF 

00 nF 



00 F 
00 hJ 



00 aF 



where q. Is the address of this routine, n is the number of intervals to fee 
examined, h is the length of an interval, and a is the address of the closed 
subroutine which replaces x in ^ by f (x) in H^ If the first parameter is 
JO kF, the roots will he printed with sign in a column to k decimal digits. 
If it is 50 rF, the roots will he stored sequentially starting at memory position 
r. If x £ be the roots found, then ,...,< x^ < * ± < » 1+1 < *•••> always. 
Control will he returned to the right side of p+3 after the n intervals hare 
been examined. 

NOTES 

If it is never desired to print the roots, the last 18 words 

of this routine (62 - 79) may he overwritten. Farther, if for other purposes 

it is desired to use Code P2 (Decimal Fraction Print), the one incorporated in 

this program may he utilized in the usual fashion by merely transferring con- 

trol to the left hand side of the 62nd word. 
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LOCATION 


ORDER 


... 
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H 2 




00 .K(H2) 











40 10F 
S5 F 




Store x~ 
Move in link 




1 


L4 6lL 
42 10L 




Plant n's address 




2 


L4 6il 
42 12L 




- Plant h's address 




3 


l4 6il 

42 8L 




Plant call for auxiliary routine 
address 




4 


42 16L 
36 7L 




Plant link address 




5 


46 59L 

T C C-7T 




Set store or print orders : 




6 


42 53L 
22 8L 








X 


46 55L 
L5 54L 


F|°om 4 






8 


42 53L 










L5 (p+3)F 


By 3 


From 6 




9 


46 18L 

46 21L 








10 


46 32L 
L5 (p+l)F 


By 1 






11 


10 20F 
01 20F 




Plant step counter 




12 


40 8F 










L5 (p+2)F 


By 2 


h to 9 




13 


40 9F 
L4 10F 




Set first interval end points 




14 


40 12F 










L5 8F 


From 53 9 


57, 60 




15 


LO 6lL, 
40 8F, 




Test for end 




16 


36 17L 









LOCATION 


ORDER 
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22 (p+3)F 


By k 


Link 


17 


L5 10F 
50 17L 


From l6 




18 


26 ( )F 
10 IF 


By 9 




19 


kO IIF 
L5 12F 






20 


40 7F 
50 20L 






21 


26 ( )F 

10 IF 


By 9 




22 


kQ 13F 
L3 13F 






23 


36 ^9L 
50 13F 






2k 


75 IIF 
36 k&L 






25 


L5 IIF 
LO 13F 


From ^-5 




26 


kO ikF 
50 10F 






27 


75 13F 








kO 16F 




Search for a root 


28 


SI F 
50 12F 






29 


7^ IIF 
LO 16F 






30 


66 ikF 
S5 F 






31 


^0 lij-F 
50 31L 






32 


26 ( )F 

10 IF 


By 10 




53 


ko 15F 




- 











LOCATION 


ORDER 
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L3 15F 






34 


36 47L , 
L5 10F 




■ 


35 


LO 12F. 
40 l6F 






56 


51 22L . 
00 2F 






37 


L2 l6F 
36 47L. 






38 


50 15F 
75 11F 






39 


32 4lL 

T C "1 1 T71 






40 


10 IF 
40 11F 






in 


22 43L 








L5 12F 


From 39 




42 


40 10F 
L5 13F 


/ 




43 


40 11F , 








L5 l4F r 


From 4l 




44 


40 12F 
L5 15F 






45 


40 13F 1 
26 25L. 






46 


80 F 




- 




00 F 




| —1 


47 


L5 14F 
22 49L 


From 34^ 


37 


48 


L5 46L 
22 49L 


From 24 




49 


L5 12F 


From 23 




| 


40 i6f 




Result in l6 




' 







LOCATION 



50 

51 
52 

53 

5^ 
55 
56 
57 
58 

59 
60 
6i 
62-79 



ORDER 



L5 TF 
40 10F 
IA 9F 
40 12F 
LT i6f 
32 53L 
22 14L 
26 ( )F 
L5 l6F 
50 5^-L 
40 ( )? 
L5 55L 
IA 28L 
46 55L 
22 l4L 
00 58L 
92 129F 
L5 l6F 
50 ( )F 
50 59L 
26 62L 
22 14L 
00 F 
00 IF 



By 6, 8 
From 53 

By 7 



From 53 



By 5 



NOTES 



?AGE 4 H~2~ 



Advance test interval 



Test result _, if < repeat 



From 50 

Store root and advance address of 

store order 



Line feed and print root 



Unit of count 

Code P2j, Decimal Fraction Print 



